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We examine the extent to which it is possible to realize the NMSSM "ideal Higgs" models es- 
poused in several papers by Gunion et al. in the context of partially universal GUT scale boundary 
conditions. To this end we use the powerful methodology of nested sampling. We pay particular 
attention to whether ideal-Higgs-like points not only pass LEP constraints but are also acceptable in 
terms of the numerous constraints now available, including those from the Tevatron and B-factory 
data, {g ~ 2)^ and the relic density Q,h^ . In general for this particular methodology and range of 
parameters chosen, very few points corresponding to said previous studies were found, and those 
that were found were at best 2a away from the preferred relic density value. Instead, there exist 
a class of points, which combine a mostly singlet-like Higgs with a mostly singlino-like neutralino 
coannihilating with the lightest stau, that are able to effectively pass all implemented constraints in 
the region 80 < nih < 100. It seems that the spin-independent direct detection cross section acts as 
a key discriminator between ideal Higgs points and the hard to detect singlino-like points. 
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I. INTRODUCTION 

As the LHC begins its first physics runs, more scrutiny can be placed on possible regions where Beyond 
the Standard Model (BSM) physics can exist. One intriguing possibility is the "ideal Higgs" scenario pointed 
out by Gunion et al. [I]-(3]. In this Next-to-Minimal Supersymmetric Model (NMSSM) scenario, the lightest 
Higgs, hi, has mass to^^ ^ 100 GeV and SM-like couplings to WW, ZZ, but decays in such a way that LEP 
limits are obeyed. In general, this is possible if the hi decays primarily to a pair of the lightest pseudoscalar 
Higgs bosons, hi — > aiOi, (resulting in a small branching ratio for hi — ^ bb). Because m^j < 2tob, the oi 
then subsequently decays to 2t, 2g, or 2c, the precise mixture depending on tan /3, with 2r being dominant 
at high tan/3 while all states are important at low tan/3. This allows these points to sidestep the LEP results 
g] on C^2b and C^l^, defined as: 

C^^ = \Cv{l)\^ X BR{hi^bb), 

C^^ = \Cvil)fxBR{hi^aiai)x[BR{ai^bb)]^, (1) 

where |Cy(l)p = gfi^zzl dhsMZZ^ where hsu denotes the Higgs boson of the Standard Model. Roughly 
speaking, LEP constraints are evaded if BRihi — > aiai) > 0.7. 
There are three main motivations for such a scenario. 

• Precision electroweak (PEW) data is most consistent with a SM-like Higgs having mass below 100 GeV. 

• There is an excess in the combined LEP data for e+e^ — > Zhb for A/^^ in the region 80— 100 GeV that 
is consistent with BR{hi — bb) < 0.3, as obtained for BR{hi — >■ aiOi) > 0.7. 
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• A mass < 100 GeV is consistent with a light superparticle spectrum for which GUT scale param- 
eters need not be fine tuned in order to obtain the correct value of mz at low scales. 

One of the primary goals of this paper is to elucidate the extent to which the emergence of ideal-Higgs- 
like scenarios depends on the extent to which the pattern of soft supersymmetry breaking is highly con- 
strained/universal. 

Indeed, a key uncertainty in both the MSSM and NMSSM is the pattern of soft supersymmetry breaking, 
as described by the scalar masses mo, gaugino masses toi/2 and trilinear couplings Aq. These presumably 
originate from physics at some high-energy scale, e.g., from some supergravity or superstring theory, and 
then evolve down to lower energy scale according to well-known renormalization-group equations. What is 
uncertain, however, is the extent to which universality applies to the scalar masses mo for different squark, 
slepton and Higgs fields, the gaugino masses mi/2 for the SU{3), SU{2) and U{1) gauginos, and the trilinear 
couplings ^0 corresponding to different Yukawa couplings. Certain types of universality are much better 
motivated than others. 

The suppression of flavour-changing neutral interactions suggests that the mo may be universal for different 
matter fields with the same quantum numbers, e.g., the different squark and slepton generations. However, 
there is no very good reason to postulate universality between, say, the spartners of left- and right-handed 
quarks, or between squarks and sleptons. In Grand Unified Theories (GUTs), there must be universality 
between fields in the same GUT multiplet, e.g., w^, d^, w_r and eji in a 10 of SU{5), and this would extend to 
all matter fields in a 16 of 5*0(10). However, there is less reason to postulate universality between these and 
the Higgs fields. Nevertheless, this extension of universality to the Higgs masses (UHM) is often assumed, 
resulting in what is commonly termed the constrained MSSM (CMSSM) or constrained NMSSM (CNMSSM). 
Alternatively, there may be non-universal Higgs masses (NUHM) in the more general MSSM and NMSSM. As 
regards the A parameters, in the MSSM or CMSSM context it is primarily the 3rd generation At that matters 
and so universality or lack thereof does not have significant phenomenological impact. However, in the 
NMSSM there are two new A parameters: one, Ax, associated with the singlet-Higgs-Higgs interaction at the 
superfield level; and a second, associated with the singlet-cubed superpotential terms. Currently, there is 
no reason for these to have the same value as Af. Indeed, the limit in which A„ and Ax are zero at the GUT 
scale is one in which the NMSSM has an extra U{1)r symmetry, independent of the GUT-scale value of At [S]. 
Although it can certainly be the case that the universalities assumed in the NUHM relaxation of the CMSSM 
and the NUHM plus Af^,Ax relaxation of the CNMSSM are still more restrictive than suggested by many 
schemes of supersymmetry breaking, comparing results incorporating such relaxations to those obtained using 
the most restrictive CMSSM or CNMSSM universalities provides a simplified framework for understanding 
what new phenomena arise when poorly-motivated restrictive boundary conditions are relaxed. In particular, 
we will see that in the NMSSM context relaxation of the CNMSSM boundary conditions to include NUHM 
and non-universality for Ax and A^ relative to At already gives rise to dramatic new physics possibilities, 
including the possibilities of ideal-Higgs-like scenarios and singlino-singlet dark matter scenarios. If nothing 
else, this shows that overly restrictive universality assumptions should be avoided in order that dramatic 
new physics scenarios are not "needlessly" excluded. 

As referred to above, previous studies in the context of the Constrained NMSSM (CNMSSM) [SlITj, in 
which all A parameters and all soft-SUSY-breaking masses-squared are unified, did not find ideal-Higgs-like 
points. In this paper, we allow the Higgs soft masses-squared to be independent of the other (unified) 
soft masses-squared and we allow the soft-SUSY-breaking A^, parameter associated with the singlet field to 
range freely, independent of the other A parameters (which are taken to be universal) . Using the advanced 
scanning techniques from SuperBayeS coupled to MultiNest [S] we can efficiently scan for interesting regions 
of parameter space and see how often points of an ideal-Higgs nature are found and look at the many 
phenomenological aspects that these points entail. In addition, the power of this approach to scan over 
all interesting parameters instead of fixing crucial ones to some canonical value can provide us with some 
insights into the full structure of the parameter space. In particular, we will not be using fine-tuning as 
the defining criteria when searching for ideal-Higgs-like points. This said, our choices of parameters are 
influenced by UrtSj and the findings therein. 
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II. THE IMPLEMENTATION OF NMSSM MODEL. 



The NMSSM superpotential contains a new superfield S which is a singlet under the SM gauge group 
SU{3)c X SU{2)l X U(1)y- (For simplicity, we use the same notation for superfields and their respective 
spin-0 component fields.) 

W = (r„ Hi u + Yd H\ d + Y^H^ V e) - e,, X S H^Hi + ^ , (2) 

where Hj = (i?^, H^^), ijj — {H^, i/"), i,j are SU{2) indices with ei2 — 1, while A and k are dimensionless 
couplings. With the addition of a scalar singlet superfield field, there will be five neutralinos and the Higgs 
content of the NMSSM is extended to include three scalar Higgses, hi, /12 and /13, and two pseudoscalars, oi 
and a2. The lightest Higgs hi plays an important role in the scenarios we consider here, and in particular 
the composition will be commented on. The state composition can be written as, 

hi = SuHu + SdHd + SsS. (3) 

The superpotential in Eq. ([2| is scale invariant, and the EW scale will only appear through the soft-SUSY- 
breaking terms in >Csoft, which in our conventions is given by 

— £soft — Q* Q + u* u + d* d + mj^ L* L + rnj^ e* e 
+ H*d Hd + ml^ H^H^ + m| S* S 
+ (a„ r„ HiQ'u + Ad Yd H'd d + A, Y, H^ U i + H.c. 

-e^jXAxSH^Hi + ^nA^ S^ + H.c, 
- ^ (Afs A3 A3 + M2 A2 A2 + Ml Xi Xi + H.c.) . (4) 

When the scalar component of S acquires a VEV, s — (S), an effective interaction —EijiiH^Hl is generated, 
with /u = As. 

In addition to terms from >Csofti the tree-level scalar Higgs potential receives the usual D and F term 
contributions: 



2 I „2 „ „2 

2/'|rr|2|c|2 , \ tj |2|c|2 , I, ijl tj] a\ , I,j2|c|4 



Vp = |A|2 {\Hdm^ + li/^n^r + \e.,HlHi\^) + l^^^l' 

-{e^jXn*WdHiS*^ +YL.c) . (5) 

Using the minimization equations we can re-express the soft breaking Higgs masses in terms of A, k, A\, 
Vd = (-ffd)i ""u = {H°) (with tan/3 = Vu/vd), and s: 

mjj^ = - A^ (s^ -I- sinV) - cos 2/3 -f As tan/3 (ks -I- Aa) , (6) 



m 



H^=- X'^ (s^ -I- cos^ /3) + ^Mf cos 2/3 -I- As cot (3 {ks + Ax) , (7) 
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TOc = - A^w^ - 2k^s^ -f Xkv'^ sin 2/3 + ^^^^ sin 2/3 - nA^s . (8) 

2s 

We are now looking at a Lagrangian that is identical in structure to our previous work fB^ but with a few 
important differences in terms of unification. In contrast to where we took CMSSM-like boundary con- 
ditions, in the present paper, we allow the Higgs mass parameters rriH^ and to freely vary, in a similar 
manner to 0. In addition A^ is no longer taken equal to the universal value, of the other A parameters. 
This freedom, specifically allowing j^^l to be small, will make it possible to obtain lighter pseudoscalar and 
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scalar Higgs masses than in ^ that are nonetheless still allowed by collider constraints. This additional free- 
dom in Higgs mass can be seen by looking at the tree level pseudoscalar mass matrix in the basis S) [TUj : 



sin 2/3 



(ks + Ax) Xv{Ax-2ks) 



\v{Ax-2k.s) X (2h + v'^ sin2l3 - SkAkS 



(9) 



After diagonalization of M^, there will be two mass eigenstatcs. The lightest state will be a mixture of the 
CP-odd doublet state Amssm that is present in the Minimal Supersymmetric Model and the new CP-odd 
component, ^5, of the complex scalar S field. We write 

ai = cos 9 aAmssm + sin6I^As, (10) 

where the entries in are such that the 11 entry is the MSSM diagonal entry. From Eqs. Q and ( 10 ) it is 
clear that having the freedom to vary is crucial if we want to control the mass of the singlet component 
independently of other parameters. Despite these changes, just as in jS] the minimisation equations, ([6|-(8]) 
will be used to fix ms, k, and s giving us a model with input parameters mi/2, ttiq, tuh^, fnHdJ Aq, A^, 
tan /3 and A, in addition to sgn(/i). In particular note that in our procedure, the value of k is an output that 
depends on A,^. It is important to emphasize that in contrast to previous studies looking at the ideal-Higgs 
region, where scanning was done at the EW scale and run up, here all of our parameters (excluding A) 
are searched over at the GUT scale and then run down. One implication is that in the studies we present 
here it is impossible to obtain the values of the gaugino masses in [T] for a given parameter point as we are 
constrained by unification considerations. That said, our choice of parameters, although constrained, leads 
to a scan that is practicable and is a useful starting point to perturb from in order to better satisfy, for 
example, fine-tuning or phenomenology. 

We also present the neutralino sector since the lightest neutralino will, by assumption, play the role of 
dark matter. The mass term in the Lagrangian is given by 

>C£L--^(*TA^x°*°+H.c., (11) 
with Ai^a given by a 5 x 5 matrix in the basis {B, W, Hu, Hd, S), 

/ Ml -Mz sinBw cos 13 Mz sinOw sin 13 \ 

M2 Mz cos ew cos p -Mz cos 9w sin (3 

-Mz sinOw cos /3 Mz cos 6w cos /3 —Xs —Xvu , (12) 

Mz sin6\Y sin 13 — M z cos 8yy sin (3 —Xs — 

\ —Xv^ —^'>}d 2ks ) 



where Mi (M2) denotes the soft mass of the bino (wino) and By/ denotes the weak mixing angle. After 
diagonalization the lightest neutralino xi (which we will denote as x from now on) can be written as: 

X = f^BB + NwW + N^Hu + NaHa + N,:S. (13) 

Finally, we note that couplings of the Higgs bosons and x to SM particle states depend upon the compo- 
sitions of the former. In particular, the coupling of h\ to WW^ ZZ is given by CyiX) = 9h{ww I QhsMWW = 
S'uSin/S + Sd cos f3 (with analogous results for /12 and /13) and the coupling of ai to bb is given by 
Ca.bb = 9a,bb/9hsMbb = COS 61^ tan/3. 

III. OUTLINE OF THE METHOD 

Following the discussion of section |lT] in this constrained version of the NMSSM the free parameters are 
given by 



9 = (mi/2,mo,TOH„,m_ff^,^o,^K,tan/3, A) 



(14) 
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CNMSSM parameters 6 


SM (nuisance) parameters tp 


10"^ < mi/2 < 4000 GeV 
10"^ < mo < 4000 GeV 
10"^ < mH„ < 4000 GeV 
10"^ < m/frf < 4000 GeV 
\Ao\< 100 GeV 
\A^\ < 10 GeV 
2 < tan /? < 20 
10^^ < A < 0.7 


160 <Mt < 190 GeV 
4 < mbinibf^ < 5 GeV 
0.10 < a,{Mzy^ < 0.13 



TABLE I: Initial ranges for our basis parameters m = {0,tp). 



Without loss of generality, one can choose A > JJ]- However, we will also (as in our previous work) fix 
sgn(/i) = +1 and then n = Xs implies s > 0. The "nuisance" parameters are treated in the same manner as 



in our previous work [6], and are shown in Eq. (15): 



= {Mt,mt{mtf''',a,{mz)^'^). (15) 
Using notation consistent with previous analyses we define our eleven dimensional basis parameter set as 

TO=(0,V), (16) 

all of which will be simultaneously scanned over. For each choice of m a number of collider and cosmological 
observables are calculated. These derived variables are denoted by ^ — {£,i,£,2, ■ ■ ■), which are then compared 
with the relevant measured data, d. 

In this study we will be using the "nested sampling" method [TT] as implemented in the MultiNest [T2] 
algorithm to efficiently explore the likelihood space. Generally speaking, for this study we will be looking at 
points of interest irrespective of statistical considerations. MultiNest provides an extremely efficient sampler 
even for likelihood functions defined over a parameter space of large dimensionality with a very complex 
structure. (See, e.g., Refs. [9} [13].) 

It should be emphasised here that although nested sampling was used to obtain the points, we will be 
drawing no statistical inferences from the results. What in effect we are doing here is using the useful 
properties of the technique (especially fast scans of high dimensionality parameter spaces and the ability to 
simultaneously scan in all parameters) to get a sample of representative points for the part of the parameter 
space we are investigating. What is presented below then, must be viewed with some observations in mind. 
Firstly, although data is implemented in the scans, the trading off of poor fits in one variable for good fits 
in another can lead to some outlandish values for key phenomenological values. In addition as the chain 
has to start somewhere, an unweighted scatter plot such as we will be showing below will also show these 
initial, poorly fitting points. To address this, for the key points we will state more precisely what the key 
phenomenological values are, and in general they correspond closely to experimental values. Also, regions 
that have clearly unacceptable experimental values will be identified where possible. 

The specific region we are investigating is defined by our range of priors which are specified in table [T] 
The above choice of priors was influenced by previous work and the preference to focus on lower values 
of the soft masses in order to explore ideal-Higgs-like scenarios. This also leads us to choose a log prior, 
defined here as being flat in logTOi/2, logmoi logm^f^, and logmn^ and flat in Aq, A and tan/3. For the 
nuisance parameters we use flat priors (although this is not important as they are directly constrained by 
measurements) and apply Gaussian likelihoods representing the experimental observations (see table |ll]), as 
before [mUTB]. 

The region specified by this set of priors is by no means a fair and even-handed exploration of the parameter 
space, but the aim here is to try and find points in a particular regime. In particular, allowing to be 
independent of Aq (unlike in our previous paper) was very important, in addition to having non-universal 
soft Higgs masses (m//„ ^ n^Ha ^o)- With this additional flexibility compared to the CNMSSM type 
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SM (nuisance) parameter 


Mean value Uncertainty 
/i o (exper.) 


Ref. 


Mt 


172.6 GeV 1.4 GeV 
4.20 GeV 0.07 GeV 
0.1176 0.002 


|21| 
I22j 



TABLE II: Experimental mean /i and standard deviation a adopted for the likelihood function for SM (nuisance) 
parameters, assumed to be described by a Gaussian distribution. 



scan we were able to find points with the desired characteristics, although as displayed below these points 
are still in the minority. 

Two alternate exploratory scans were implemented to better understand this region, the first being to allow 
a much more generous range in the parameters, and the second was to do a scan with similar constrained 
priors but with the unification conditions = ''^Hd — "^o ^iid = Aq enforced, as in the so-called 

CNMSSM studied in The objective was to see if it was the focusing into a small region or the extra 
freedoms in the Higgs sector that lead to us finding points of interest. What we found is that both the 
freedom and the focus (especially in A^) seem to be important. 

We compute our mass spectra and observable quantities using the publicly available NMSSMTools (version 
2.1.1) that includes NMSPEC with a link to Micromegas [T7]; for details see Ref. [TH]. We hst the observables 
that the current version of NMSPEC is applying to points found in the scan in Table III The relic density 
ri/i^ of the lightest neutralino is computed with the help of Micromegas, which is also linked to NMSPEC. 
We further use the same code to compute the cross section for direct detection of dark matter via its 
elastic scattering with targets in underground detectors but do not include it in the likelihood due to large 
astrophysical uncertainties. The likelihoods for the measured observables are taken as Gaussian about their 
mean values, /i as tabulated in Table |III[ the Gaussian widths being determined by the experimental and 
theoretical errors, cr and r, respectively (see the detailed explanation in Refs. [131 [H]). In the case where only 
an experimental limit is available, this is given, along with the theoretical error. The smearing out of bounds 
and combination of experimental and theoretical errors is handled in an identical manner to Refs. |14| , 115 ) . 
with the notable exception of the Higgs mass and LEP limits on sparticle masses, which are calculated as a 
step function with values of the cross section times branching ratio (in the case of the Higgs) or mass that 
are within two standard deviations of the experimental limit being accepted. Finally, any points that fail to 
provide radiative EWSB, give us tachyons or the LSP other than the neutralino are rejected. 

The above discussion does not yet include the constraints from the recent analysis by ALEPH for the 
e+e" — > Zh, h aa At channel [inj nor the constraints from BaBar data on T{3S) — 7a with 
a — > T+r~ [20]. These will be considered ex-post-facto. In this way, we can see what the impact of these 
latter constraints is upon a less biased sample of otherwise acceptable points in parameter space. 



IV. RESULTS 



In this section we present our numerical results from global scans, mostly in the form of scatter plots for 
some of the most interesting combinations of observables. Most of the figures displayed are scatter plots with 
three distinct populations visible on them. There are a large number of grey points in the figures that follow; 
these so-termed full scan points come from all the points obtained in the scan and are thinned by a factor 
of one hundred for clarity. This will also eliminate many of the outlying initial points talked about earlier. 
This population of points can be thought of as a representation of the general structure of the parameter 
space, the backdrop against which we specify points of interest. 

To better display the points with BR (oi bb) =0 and BR{hi -> aiai) > 0.5 we have marked them 
as triangles in the relevant figures. Such points will be termed Type I points. In addition, there are 
points with BRihi — >■ bb) < 0.5 by virtue of substantial BR{hi aioi) but for which > 2mB and 
BR (ai — > bb) 7^ 0. These are shown by squares in the figures. We call these points Type II points. In 
fact, there are two subclasses of Type II points. Type HA points are such that the light Higgs is mainly 
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0bs6rv3,blG 


Nlean value 


Uncertainties 








(T ( f^~^Tnf^r 






S(n — 2) X in^" 


29.5 


8.8 


1 


I23l 


Xs-t) X 10* 


3.55 


0.26 


0.21 


|23| 


ru) X 10* 


1.32 


0.49 


0.38 


m 




0.1099 


0.0062 


0.1 n/i^ 


|25| 




Limit (95% CL) 




r (theor.) 


Ref. 


BRiBs M^M") 


< 5.8 X 10"** 




14% 


\H 




As implemented in NMSSMTools. 






|18| 


sparticle masses 


As implemented in NMSSMTools. 






|18| 


'-'eff 


As implemented in NMSSMTools. 






|18| 


^46 


As implemented in NMSSMTools. 






|18| 



TABLE 111: Summary of the observables used in the analysis. For more details on how these are implemented, see 

mi- 



doublet. Type IIA points pass constraints on C^g and C*^ considered individually, but may struggle to 
pass the overall constraint implicit in LEP data on hZ production with h — 6's, where h' s represents any 
final state containing one or more h quarks. This overall constraint becomes important for below about 
110 GeV. Because of the need to consider this overall constraint, one cannot be certain of whether or not 
the Type IIA points with m/j^ < 110 — 112 GeV should be eliminated without submitting them to the LEP 
collaborations for full analysis. Thus, we will depict them in the figures. However, in general the Type I 
points where < 2tob are more interesting. Type IIB points are ones with BR{hi — >■ bb) < 0.5, large 
BRihi — 7> aifli) and > 2mB but for which the light Higgs is mainly singlet. These easily pass the LEP 
constraints and have other interesting properties that we shall elucidate later. 

There is a third class of points. Type III, that are worth singling out. These have BR{hi — > hb) > 0.5 but 
pass our primary selection criteria. As for Type II points, there are two subclasses: Type IIIA for which the 
hi is mainly doublet and is small; and Type IIIB for which hi is mainly singlet, S'^ ~ 1. The latter easily 
pass LEP limits on C^g and C^g, by virtue of hi being mainly singlet in composition, implying very small 
ZZhi coupling. In common with Type IIB, the Higgs with SM-like coupling to WW, ZZ is instead the 
which must have mass and/or decays that allow it to obey all of the bounds imposed by LEP. Type IIIA 
points pass LEP limits by virtue of m^.^ > 114 GeV, and are therefore not "special" in any way. Type IIIB 
points have € [80 — 100] GeV and are particularly interesting in that they have a high probability of 
providing the correct value of flh^ . In fact, in what follows we will include in the definition of Type III points 
the requirement that they yield the correct value of flh^ within 2a. This will mean that Type HI points will 
have the lowest overall of all the points in our scans. However, they cannot explain the LEP excess in 
this range since, as stated before, the hi is for the most part decoupled and it is the which plays the 
role of the SM-like Higgs, and it has mass above 110 GeV. Finally, we also include in our final definition 
of Type HI points the requirement that they satisfy the "ex-post-facto" constraints from ALEPH (on the 
hi flifli 4r channel) and from BaBar (from T{3S) jai decays) mentioned earlier. In short, our final 
definition of Type HI points is such that they are very consistent with all available experimental constraints. 



We explicitly state our selection criteria for the various points in Table IV 



A later table will provide a fuller list of the properties of the different classes of points. 

Given the motivation for searching in those regions of parameter space that might yield ideal-Higgs sce- 
narios, we first look at the hi branching ratios and C^g. The first result is shown in Fig. Ill where we plot 
C^g as a function of m/i^ for our study, points from previous work by Gunion et al. and the experimental 
limits from LEP. This plot already shows that the only points in our scans that reproduce exactly the desired 
qualities {BR (/ii 66) ^ with BR (ai hb) — 0) while passing all LEP constraints are precisely those 
with a Higgs mass in the range 80 GeV < ruh^ < 100 GeV. 

In Fig. [Tl one also clearly sees the Type IIA points with > 108 GeV and the Type IIB points with 
"^/ii % 90 TleV. As stated above, those of the Type IIA points that have rn/n < 114 GeV, especially those 
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Point Type 


BR {hi bb) 


BR 


(ai 


bb) 


5? 


Type I 


< 0.5 




= 




~ 


Type IIA 


< 0.5 




/o 




~ 


Type IIB 


< 0.5 




/o 




~ 1 


Type IIIA 


> 0.5 




/o 




~ 


Type IIIB 


> 0.5 




/o 




~ 1 



TABLE IV: Definition of various points shown in figures. In the case of Type III points, we have also required that 
flh^ be within ±2cr of the observed value and that they obey the ALEPH constraints on hi — >■ aiai — )■ 4r and the 
BaBar constraints on T(3S) — > 701. 
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FIG. 1: A plot of the LEP limit from g] in the C^g vs TUhi plane, superimposed on points obtained by Gunion et 
aim [1] and so called "Type I" points {rriai < 2mfl) and"Type 11" points {niai > 2mB) from this study. 



with m/ij < 110 GeV, might be ruled out by a combined Z + b's LEP analysis, even if not ruled out by the 
C^g and separate limits. 

We display in Fig. [2] an expanded look at this region. In the top left plot we put into context the correlation 
between nihi ^^'^ for Type I and Type II points by showing "background" points from the full scan. 
There, we see that there are many background points with C^g of order 0.7 to 0.8 but with large enough 
m/ij to escape LEP limits, as characteristic of Type IIIA points, as well as background points with very low 
m/ij and C^g, which include Type IIIB points. Other plots in this figure show that we obtain the points 
with small C^g in two different ways. The first means is to suppress the branching ratio, BR (/ii -> bb), as 
is the case with most Type I and Type IIA points. The second is to suppress the squared-coupling |Cy(l)p, 
as happens if the hi is sufficiently singlet in composition, which is the case in particular for Type IIB 
points. Indeed there can be a lot of points in a similar m/j^ region to that identified by Gunion et al. {i.e. 
80 GeV < m/ij < 100 GeV) that escape LEP limits via suppression of |Cv/(l)P rather than via suppression 
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FIG. 2: Plots of various phenomenological quantities — see Eq. ([TJ — as a function of . Note that all plots and 
all figures to come use the color red for the Type II points, as opposed to the coloring in Fig. [T] where the color blue 
was used for Type II points. 



of BR{hi — > bb). Small |Cv(l)p implies that the hi cannot act as a "ideal" Higgs defined as having SM-like 
WW, ZZ coupHng but mass < 105 GeV. 

The bottom right figure in Fig. [2 shows C^^ against m^j^ . This can be compared with a similar figure in 
[T]. In general it is clear that the Type I points all have C^g = (since < 2m b) while for Type IIA 
points C^g is quite significant and, for those points with < 114 GeV, is not far below the LEP limit. As 
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noted earlier, since the limits on C^^ and C^g are being applied individually and not in combination it could 
well be that the Type IIA points with m/^ < 114 GeV, especially those with m/^ < 110 GeV, are in fact in 
contradiction with LEP. But, without a full LEP analysis, it can be instructive to leave them in with this 
caveat in mind. One thing to notice in general is that the ideal-Higgs-like Type I and the Type IIA points 
with m/ij < 110 GeV and fairly large |C\/(1)| ~ 1 are very rare, even in the context of a scan looking for 
these regions; this could be an artifact of our scanning technique or it could be that these points are truly 
hard to find given the criteria and high-scale boundary conditions we have used. 
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FIG. 3: Plots of TTiai as a function of rrihi and of C^g as a function of . 



Since is so crucial to whether or not a given point is ruled out by LEP data, it is useful to understand 
how mill s-iid C^g correlate with ■ In the left plot of Fig. [s] we display all the different types of points 
in the TO^j-mai plane. Note again the Type IIA points with rua^ > 2mB and 109 GeV < m/^ < 114 GeV 
that escape LEP limits on despite having large BR{hi — > aiai) and large BR{ai — > 66). These also 
appeared in Fig. [2] 

To better demonstrate the interplay between the various branching ratios needed to evade LEP constraints, 
some of the crucial ones are depicted in the same convention in Fig. |4j One can see that the green Type I 
points are clearly isolated, with the key discriminator from Type II points being BR (ai — )■ 66) . In Fig. |4j 
some key differences between Type IIA and Type IIB points are apparent, the most notable being the very 
small BR{hi — ^ aiai) for Type IIB (singlet hi) points. Note also that all Type III points have very small 
BR(hi — )■ oiai). 

We will shortly discuss whether or not the Type I points escape the latest ALEPH limits on hi — > aiCi 
with fli — > T'^T^. Such escape is possible when tan/3 is small, since at small tan/3 one predicts that 
BR{ai — > T+T^) is significantly suppressed due to substantial branching ratios for ai to cc, ss and gg and 
the resulting final states in hi — >■ aiOi are less strongly constrained than the hi — >■ aiai — > 4r final state. 
This was discussed in [57] . 

In Fig. [5] we show the square of the singlet component of the /ii, and the square of the singlino component 
of the lightest neutralino, the Xi ^ functions of ■ These figures illustrate a number of things. First note 
the large number of points with iV^ ~ 1 and ^ 1, the latter implying that |Cy(l)p is greatly suppressed. 
Included in this set of points are the Type IIB points with nihi < GeV as well as the Type IIIB points 
with low m/ij, large BR(hi — ^ aiai) and non-zero BR{ai — >■ 66) — all these points escape LEP limits since 
the hi is very singlet-like. Second, we observe that the S'^ plot is closely related to the top right panel in 
Fig. [2j As noted in the discussion of the latter figure, it is the Type I points and Type IIA points with 
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FIG. 5: We plot the square of the singlet component of the hi and the square of the singlino component of x as 
functions of nihi- 



rahi > 80 GeV for which the hi is highly doublet-like whereas the Type IIB points with nih^ < 80 GeV 
have a singlet-like hi. And, finally, there is the large collection of points with m}-^^ > 114 GeV that are 
unconstrained by LEP data and typically are very doublet-like. 

The right hand plot of Fig. [sjwill be of more use below when we consider dark matter, but it does show an 
interesting correlation, namely that the Type I points and the Type IIA points with m/n > 80 GeV stand 
out by having a x that is bino-like instead of singlino-like, the latter being more typical of the majority of 
points found in our scans, including Type IIB points. We further note that the Type IIIB points that have 
nihi < 110 GeV and ~ 1 also have N"^ ^ 1. In contrast, the Type IIIA points which have > 114 GeV 
and small S"^ (with being large instead) can have either large iV^ or large {i.e. singlet-like x or bino-like 
X)- 
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1 


2.35 


1.78 X 10"^° 
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BR (hi — )• aiQi) 


BR {hi bb) 
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I? lALEPH 
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1 


0.341 
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0.0943 


0.833 
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.2902 


0.615 ! 


2 


0.199 


0.0197 
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0.844 


0.128 
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0.0042 


.49809 


0.655 ! 
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0.664 


0.0136 


0.890 


0.0893 


0.771 


0.0047 


.75916 


0.529 


4 
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0.0133 
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0.916 
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0.786 


0.0034 
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TABLE V: Displayed are some values of interest for the Type I points found in our scans. In the upper table we show 
the base parameters that give us our population of interesting (Type I) points. The final column, denotes G, defined 
in Eq. ( 18 1, a measure of the fine-tuning needed to obtain the (low) value of ruax ■ the middle table are some of the 
phenomenological values for the points of interest. Notice the likelihood (to be precise the — 21og(likelihood) = x^) 
in general is large refiecting a poor fit, and this is largely being driven by poor fits to Q,h? [13]. The bottom 
table shows some of the key branching ratios of interest for Type I points, and compares the ALEPH limits on 

) X [BR{ai — >■ T^T ] with the predicted values. Points appended with an exclamation 
mark are excluded by the ALEPH analysis. However, as discussed later, by adjusting A^, by a very small amount 
they can be brought into agreement with the ALEPH limits without affecting any other phenomenology. 



Full details regarding Type I points appear in Table [V] The upper table shows the input parameter values 
for each of the Type I points. The corresponding "light-ai" fine-tuning measure, G, defined by: 

G = Mtn{[Maxi\FAjAFAM\FA, + FaJ} , (17) 

where, 

is also shown. One can see some common threads for all of the Type I points. Perhaps the most intriguing 
is the need for A to be quite large and away from the decoupling limit. This suggests that these points are in 
some sense specific to the NMSSM and are unlikely to be found in similar parametrizations of the MSSM. It 
is also nice to see that despite the light-ai fine-tuning measure G not being used in the scans, the resultant 
values for Type I points are not wholly unreasonable. Finally, we note that the values of tan /3 for which 
Type I points were found are relatively low. In the less constrained scans of parameter space performed in 
[TH3] Type I (ideal-Higgs) points were found at large tan /3 as well. 
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In Fig. IgI we show the values of toq and vn^i^ (at the GUT scale) for the various different types of points. 
We see that many of the points of Type II and Type III have quite small values of TOq and that most Type I 
points have quite modest mo. As regards it is typically of order 250 — 300 GeV for Type I points but 

ranges from ^ 250 GeV up to 2 TeV for Type II and Type III points. Regardless, the resulting gluino and 
(non-stop) squark masses are always at least as large as 650 GeV and often significantly larger. Such values 
are above the limits currently being set by LHC data, which limits are typically of order 500 — 600 GeV 
(assuming universality for the gaugino masses and for the non-Higgs scalar masses at the GUT scale). Of 
course, the LHC will probe gluino and squark masses of order 1 TeV after another year or two of running. 
In common with other models employing universality at the GUT scale, the parameter points typical of our 
study will then start to be ruled out. 

The middle table in Table [V] gives some corresponding experimental values for the Type I points. It is 
interesting to see that the points provide phenomenologically viable results for BR{B — )■ ^^7) and BR{Bg — 
fi^fi^). However, the lower 2a boundary for the observed 6{g — 2)^ is ~ 10 x 10^^'^ and only Type I point 2 
barely predicts this high a value, the other Type I points predicting values in the range (1.78 — 3.90) x 10~^°. 
The relic density is equally problematical, with the best value barely getting to within 2a of the WMAP 
value. The likelihood is dominated by this contribution as in general the relic density is the strongest 
constraint on the parameter space [13]. Hence, the points with best likelihood correspond to ilh^ closest to 
its experimental value. As discussed shortly, the Type IIIB points {i.e. points nih-^ < 104 GeV that escape 
LEP limits by virtue of the hi being mainly singlet) quite readily achieve an fi/i^ near the WMAP value; as 
a result, the for Type IIIB points ranges from a low of ~ 1.9 to a high of ^ 6, vs. the best value of ^ 26 
found for Type I points in our scans. Type IIIA points have values only slightly larger than Type IIIB 
points and definitely below 26, as consistent with our requirement that Type HI points be consistent with 
the observed il/i^ within zL2a. Note that small can be achieved within the other defining characteristics 
for Type IIIA and Type IIIB points because small fi/i^ is possible despite the singlino or bino nature of the 
X by virtue of near mass degeneracy of the x and ti . 
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The bottom table of Table |V] gives the values of the coupling C^^f,^ = tan (3 cos 9 a in comparison to the 
maximum absolute value allowed by BaBar data in the T35 — >■ 7r + t~ channel. Here, cos 9a is the doublet 
component of ai as defined by ai — cos 9 aumssm + sm9Aas- This bottom table also gives the value of 

= |Cy(l)p X BR{hi oiai) x [BR{ai — > r+r^]^ in comparison to the upper limit for each point from 
the recent ALEPH analysis. We observe that the Type I points have no problem obeying the limits from 
BaBar but that the ALEPH limits are very problematical for three out of five of the Type I points. However, 
we show below that a very small change in A^, will bring these points into agreement with the ALEPH limits 
without affecting any other phenomenology. We also wish to note that the ALEPH constraints are much 
stronger than what was expected on the basis of Monte Carlo and so, in our opinion, some relaxation of the 
ALEPH bounds could be considered. If ^ la relaxation is allowed, then all our Type I points survive "as is" . 
We also wish to note that ^- x BR{ai — > jjT) roughly determines the ability to detect gg ^ ai ^ M"*^/^" 

at hadron colliders [53]. Very roughly, in the niai < 2mB region where BR{ai — > ji^ ^i^) ^ (0.003 ~ 0.005) 



detection will only be "easy" if IC^^bjl > 1. 
with point 2 providing the largest value of ^ 



Unfortunately, for our Type I points, 
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FIG. 7: We plot as a function of its singlet component, | cosOa\ and ^j^BR{ai — > fi'^ /i ) as a function of 

A more global picture of | cos 6*^1 and ^j^BR{ai — ^ M^A^ ) is provided by Fig. jTj In the left plot, we 
give TTiaj as a function of the magnitude of its singlet component, | cos0^|, for all points in order to show 
more generally how singlet the ai is for the different classes of points. We observe that the ai is extremely 
singlet for the vast bulk of points, including the Type IIB and Type IIIB points. This correlates with the 
fact that the hi is mainly singlet for these same two types of points. 

In the right plot of Fig. jTjwe show ^-BR{ai as a function of to indicate which points 

have a reasonable probability that the production/decay channel gg ^ ai ^ A'^M^ could be detected. Points 
for which this product is > 0.001 would have a viable signal at the LHC for accumulated luminosities of 

order 10 fb^^ 28 (more being required if is in the region of the T resonances). We see that none of 
our points are even close to allowing such detection. 

Of course, the size of C'^^^f) derives both from cos9a and tan/3. Thus, it is perhaps useful to display tan/3 
as a function of and m/^ . This is done in Fig. Is] We see that only the Type I and the Type HA points 
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FIG. 9: We plot our points in the - ruhi ^'^d show |Cv(2)p as a function of 



with mill ?i 80 GeV (for which the ai and hi have at least a modest doublet component) are forced into the 
low tan/3 region. In contrast, Type IIB and Type IIIB points all have tan /? > 8. Type IIIA points can have 
any tan/3 above ~ 3. In the left plot of Fig. [9] we show how all the points are distributed in the - m/j^ 
plane. This plot shows very clearly two branches for all the points that are neither Type I nor Type IIA 
points. The vertical branch corresponds to Type IIIA points where m/^ > 114 GeV (thereby escaping LEP 
limits) with |Cv(l)P ~ 1 (see Fig. [2|. The horizontal branch encompasses the Type IIB and Type IIIB 
points for which the hi is very singlet and it is instead the /12 that is SM-like with ^114 GeV and 
|Cy(2)p ~ 1, as displayed in the right plot of Fig. [o] Finally, we remark that BR{h2 aifli) is very small 
for all points — this is perfectly OK since LEP limits for the /i2 are evaded either because it is very singlet 
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or because > 114 GeV — the extra /i2 — >■ aiCi decay channel is not needed. 
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Let us now return to the ALEPH Umits on ^ \Cv{l)\'^BR{hi -J> aiai)[BR{ai r+r")]^. In Fig. 
we show specifically the effect of varying on One can see that there exist two "discontinuities" as 
Ak, is varied. These arise because TOqi decreases with increasing . The first abrupt change occurs at the 
upper range of plotted when niai passes below 2mr, and goes to zero because BR{ai — T+r~) = 0. 
The 2nd abrupt change arises as A^ is decreased and niai becomes larger than 2mB- In this region, the 
dominant decay channel for the lightest Higgs is hi — > aiai — ?> 46 and BRihi — > aioi — > 4t) is very small. 
As decreases just a little bit more, the model point in question will start to exceed LEP bounds on the 

in the case of point 5 the LEP bound on Cfj^ 
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hi — > 46 final state. As seen in the left hand plot of Fig. 

comes into play quite quickly as A^ is decreased. In comparison, the right hand plot shows that in the case 
of point 1 there is a larger range of A„ for which the LEP bound on Cfj^ is satisfied and = so that 
the ALEPH bound is automatically satisfied. Of course, the nominal value of A^ from Table V for point 5 
(the left-hand plot) is such that the ALEPH bound is satisfied and no adjustment of A^ is required. We 
presented the plot to show how sensitive the LEP phenomenology is to A^,- In this case, the value of can 
be changed somewhat from its nominal Table V value without immediately encountering a problem with 
either the LEP bound on the 46 final state or the ALEPH bound on In contrast, the nominal value of 
for point 1 in Table V gives a value for that is considerably too large in comparison to the < 0.29 
ALEPH limit. In this case, we must lower A^ in order to satisfy the ALEPH bound. Roughly, any A^ above 
5.58 GeV, the value at which the LEP 46 bound (dashed line) becomes relevant, but below about 5.67 GeV 
(vs. the nominal value of 5.69 GeV) would be allowed. 

The above discussion illustrates that acceptable can be obtained for the type I points 1, 2 and 4 by 
very small shifts in A^ that do not any way affect the remainder of the phenomenology of these points. Of 
course, it must be acknowledged that there is a certain level of fine-tuning of A„ involved in getting nia^ 
into an allowed range. This is, in fact, already reflected in the somewhat large G values of points 1 and 2 
of Table V, these points being ones where the nominal is substantially above the ALEPH limit. Point 4 
has a much more modest G value and correspondingly a broader range of would allow it to satisfy the 
ALEPH limit that is only slightly below the value of predicted by the nominal A„ value tabulated in 
Table V. 

As state d ea rlier, ^Ih^ plays a pivotal role in determining the likelihood of a given point in parameter 
space. Fig. 11 shows us that the range of values for Q,h^ is huge, with many points having relic densities 
that are too large by orders of magnitude. In this context, the fact that Type I points tend to achieve the 
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right order of magnitude seems quite remarkable. That said, despite the rehc density constraint pushing our 
scanning quite strongly towards the WMAP value, we did not find Type I points with a relic density that 
is less than two sigma away from the observed value. Similar remarks apply to the Type IIA points. In 
contrast. Type IIB points (for which the x is very singlino-like) have much too large Q,h^ as a result of too 
small an annihilation cross section. 

As noted already. Type III points were defined by requiring not only BR{hi — >■ bb) > 0.5 and iriai > Iuib 
(so that BR{ai — >■ bb) ^ 0) but also by demanding that Vth^ is within ±2cr of the observed value. Thus, it is 
mainly the Type III points that populate the band in the right-hand expanded plot of Fig. |11[ It turns out 
that for all the Type III points the dominant process responsible for getting correct Vlh? is coannihilation of 
the X with ti — they are quite closely degenerate in mass for all the Type IIIA and Type IIIB points. In the 
case of points where the x is very singlino-like, which comprises all Type IIIB points and a sizable fraction 
of Type IIIA points, the mass difference between x ^-nd ti is at most about 4 GeV and the common mass is 
typically of order 120 GeV. The x h&s just enough gaugino and higgsino components (of order 10~^ at the 
probability level) to couple effectively and coannihilate with the ti to get the right relic density. In the case 
of the Type IIIA points, for which the x is mainly bino-like, the common mass is most often > 250 GeV (but 
not always) and for such points the mass difference is more typically of order 10 GeV. The very smallest x^ 
values are achieved for the Type IIIB points for which both the hi and the x ^'^e mainly singlet and singlino, 
respectively, the "singlet-singlino" (or SS) scenario. The ease with which such low points were found in 
our scans suggest that the SS scenario for dark matter should be taken quite seriously as possibly being the 
correct paradigm for dark matter. 
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One question is whether small parameter changes for the Type I points could bring the predicted flh^ into 
close agreement with observation. After all, our scans might just be slightly missing Type I points with the 
right relic density. To examine this, an attempt to vary the gaugino masses independently at the GUT scale 
was made to see if this could lead to the right amount of dark matter. In general it is not difficult to get to 
ilh^ ^0.1, but in all the cases studied, it comes at the price of exceeding LEP limits on several channels, 
including Higgstrahlung processes such as Zhi production with hi — bb. 

In order to do this perturbation, we took the most promising Type I point {i.e. the one with the best relic 
density value, point 5 in table [v| ) and perturbed M2 away from universality at the GUT scale. The results 



are shown in Figs. 12 and [13] What seems to be happening here is that as AI2 changes in this region, rua^ 
also changes dramatically, so much so that perturbing M2 by only a few GeV gives us a point that is ruled 
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out by the h ^ aa ^ Ab's limit from LEP. As M2 gets bigger, eventually we get to the crucial point where 
BR (oi — >■ bb) is kinematically suppressed as ruai is sufficiently light. 
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FIG. 12: The effect of varying the M2 parameter of point 5 in table M The left-hand plot gives the values for C^ff 
as defined in Eq.Q with red denoting the experimental limit and black the NMSSM predicted value as a function 
of M2. The black curve terminates as A/2 increases when m^-^ becomes negative. The sharp rise in the experimental 
limit occurs as rriaj approaches and then falls below 2ms, the point at which the ai — > 2b decay mode becomes 
kinematically forbidden. On the right, we show the rather dramatic change of with M2. In contrast, rrihi 
remains roughly at 102 GeV over this range of M2. 



The curves in Fig.[T2]do not extend below M2 ^ 367 GeV or so, since, as seen in Fig. [13] in this region one 
exceeds the LEP limits on C^g. In addition, the sensitivity of niai means that soon after getting to a point 
where LEP limits are respected, (the exact point being unclear due to the resolution of the exploration done) 
rriai is driven tachyonic. Looking at this admittedly very specific case, one can perhaps begin to understand 
why we are obtaining so few points that are Type I and anywhere near the right relic density, given that rua-i 
is very sensitive to changes and is the crucial element here. This is by no means a blanket statement about 
points being ruled out, merely an observation for this specific situation. For further study, a scan with all 
the gaugino masses disunified might be useful to shed light on this. 

Let us now turn to predictions for direct detection of the ncutralino via scattering on nucleons. In Fig. [14] 
we show the spin-independent cross section as a function of both and , and see something not entirely 
unexpected: the large swathe of singlino points, which comprise Type IIB, Type IIIB and a sizable fraction 
of Type IIIA points, will be nearly undetectable in any upcoming direct detection experiment, having cross 
sections of at most 10^^^ pb or so. On the other hand, about half the Type IIIA points have > 10^^° pb, 
with some having (7^^ ^ 10^^ pb (and good VLh^). And, we should again note that all Type IIIA points 
have rrihi > 114 GeV. In contrast, the Type I and Type IIA points have nihi < 114 GeV and often 
TO/ij ^ 100 GeV, i.e. in the region of interest for explaining the LEP excess near 100 GeV. As apparent 
from Fig. [14] these points for which the hi is doublet-like seem to have a large, or at least measurable, direct 
detection cross section. As a result, if a largish value of (Xp ^ is eventually observed, the value of ruhi could be 
used to distinguish between the Type IIIA and Type I/IIA regions of parameter space. Note that a priori 
there is no reason for Type I and Type IIA points to have a direct detection cross section that is both high 
enough to be tested by future experiments and low enough to avoid current constraints, and as such it is 
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0.35 




&0 300 350 400 

(GeV) 

FIG. 13: The effect of varying the M2 parameter on point 5 in table |V] Here we show the change in the rehc density, 
flh^. It is possible to bring DM in line with experiment, but in doing so one violates LEP limits on several channels, 
below M2 ~ 360 GeV in this case, which is the reason that the curves in the previous figures don't go below this 
value. 
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FIG. 14: The dark matter spin-independent cross section a'p as a function of, on the left, and, on the right, 
mhi- Included on the left hand side is an illustrative limit from the Xenon-100 direct detection experiment |29| . 



interesting to note that many such points do appear. 

On another note, by comparing between the two figures in Fig. |14[ we can again detect the correlation 
between singlet-like Higgses and singlino-like neutralinos that was apparent in Fig. [5] for values of the Higgs 
mass in the range m/i^ < 90 GeV. As already apparent in Fig. [5j for 114 > m/i^ > 90 GeV the situation 
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my 



is more complicated and there is no clear correlation between the singlet nature of the hi and the singlino 
nature of the x- 

It is interesting to comment on results for (7^ ^ for points with very low < 15 GeV. We see in Fig. 
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that 



the Type IIB points are those that populate the very low region. However, these points all have a smglino- 
like X and, correspondingly, the largest cross section is of order a^^ ~ lO^^'^pb, i.e. far below the region 
that is needed to explain the possible CoGeNT [3D] and DAMA [3T] excesses in the 6 GeV < niy. < 9 GeV 
region for which > 10~^ pb is required [32]. In [33], it is shown that if one ignores GUT-scale unification 
then values of a^^ within a factor of 10 of the above range are possible while still having a Higgs with 
SM-like WW, ZZ couplings that is sufficiently light to achieve "ideal-like" agreement with precision data. 
From the study presented in the present paper, it seems that such large values of cannot be achieved in 
the context of the relaxed-CNMSSM boundary conditions employed here. In particular, very low values 
for a bino-like x (which allowed for the largest a^^ values in [33]) require small Mi values, a region that is 
quite inaccessible in the strict CNMSSM context. 

Before concluding, we provide a tabular summary of the most important characteristics of the five different 
classes of points that we have particularly focused on: Type I, Type IIA, Type IIB, Type IIIA and Type IIIB. 



V. CONCLUSION AND SUMMARY 



In the version of the NMSSM studied in this paper, we have relaxed the unification condition of and 
of the Higgs soft masses mn^ and w//^ (at the GUT scale) with respect to the CNMSSM in order to explore 
the extent to which the absence of the ideal-Higgs-like scenarios in the CNMSSM scans depended upon these 
particular (rather unmotivated) universality assumptions. Allowing for non-universal and non-universal 
Higgs soft masses, five parameter space points corresponding to the so called "ideal Higgs" scenario were 
indeed found, although these were far outnumbered by other points. The phenomenology of the ideal-Higgs 
points was studied, and in the context of this particular scan these points were seemingly acceptable in terms 
of flavour observables like BR{B Xg^) and BR{Bs ii^ fjT) (with {g — 2)^, in common with every other 
scan we have done, struggling to fit the observed 3ct difference relative to the SM). However, only two of the 
five ideal-Higgs-like points we found are strictly consistent with the latest ALEPH limits on the Higgs to 
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TABLE VI: Important characteristics of the different classes of points. The top entries in the Table list the require- 
ments imposed in defining the classes. The bottom entries list the resulting properties of points in the different 
classes. Three additional requirements are imposed for Type III points (only): (i) that ilh^ be within 2a of the 
observed value (roughly log[f2/i^] ~ — 1); (ii) ALEPH limits on hi —¥ aiai 4r are satisfied; and (iii) BaBar limits 

on T(35') — >■ 7ai are satisfied. 





Type I 


Type IIA 


Type IIB 


Type IIIA 


Type IIIB 


BR{hi bb) 


< 0.3 


< 0.5 


< 0.5 


> 0.5 


> 0.5 


BR(ai bb) 







^0 




/o 




~ 


~ 


~ 1 


~ 


~ 1 


|C^v^(l)l 


~ 1 


~ 1 


< 1 


~ 1 


< 1 


rrihi 


e [95, 108] 


> 75 


< 90, mostly < 20 


> 114 


< 110 




< 2ms 


G [2ms, 50] 


G [2ms, 20] 


> 2ms 


G [2ms, 40] 


BR{hi — >■ oioi) 


> 0.7 


> 0.6 


- 


< 0.1 


~ 


n! 


< 1 


< 1 


~ 1 


< 1 or ~ 1 


~ 1 


Nl 


~ 1 


~ 1 


< 1 


~ 1 or < 1 


< 1 


tan B 


< 3.5 


< 6 


> 8 


G [2,20] 


> 8 


cos Oa 


e [0.005, 0.02] 


G [0.003, 0.02] 


- 


G [0.002,0.017] 


~ 


m-x 


G [80, 200] 


G [70, 300] 


< 15 


G [113,400] 


G [91,110] 




G [-0.65,0.4] 


G [-1.9,1.8] 


G [6.3, 7.3] 


~ -1 


~ -1 


log[af (pb)] 


G [-8.7,-8.2] 


G [-10,-7] 


G [-14.5,-9.8] 


G [-10,-7] 


< -12 




e [26, 70] 


> 10 


> 10 


e [2.7, 26] 


e [1.9,6.0] 



four tau mode, with a third being very close to consistency. However, we have shown that by changing 
by a very small amount compared to the nominal value found in the scan (which did not use the ALEPH 
limit on the four tau final state as an input to the chi-squared employed) will allow consistency with the 
ALEPH limit without altering any other phenomenology. 

As regards the relic density, the Q.h? values of the Type I points were not within the two sigma range of 
the observations, but four out of the five were within a factor of 2 or 3 of Q.h? ~ 0.1. Perhaps not too much 
should be read into this as the relative scarcity of these points in our rclaxed-CNMSSM scan could mean our 
understanding of these parameter points is incomplete. Correct ^h? can be achieved for the Type I points 
by varying the most relevant gaugino mass parameter (M2) slightly. However, we found that for M2 such 
that the relic density was correct one or more of the LEP Higgs limits was not satisfied. 

There was a another very interesting class of points, denoted Type HI, that appeared in our relaxed- 
CNMSSM scan. Type III points are, first of all, characterized by BR{hi bb) > 0.5 (i.e. the normal SM 
decay is dominant) and by niai > 2mB (ie. ai — >■ bb is dominant). Points satisfying this criterion are, as 
one might expect, very numerous. Further, we found that it was very easy to find points satisfying the 
above criteria that gave flh^ values in close agreement (±2(7) with the observed value (something that we 
included in our final definition of Type III points, along with the requirement that they obey the ALEPH 
limits on hi aiai At decays and the BaBar limits on T(35) — )• 701 decays). Within the Type III 
class, as finally defined, the very best predictions for Vth"^ were obtained for cases in which the lightest 
Higgs is very singlet-like with 85 < rrih^ < 110 GeV and the lightest ncutralino is very singlino-like with 
91 < < 110 GeV. Such scenarios are dubbed SS scenarios. For all such SS scenarios, sufficiently small 
in agreement with experiment is achieved via X coannihilation. We believe that one should take 
these SS scenarios seriously. It will then be the second lightest Higgs boson (predicted to have mass nih^ 
close to 114 GcV) that has SM-like couplings to WW,ZZ and its decays will also be SM-like. Unfortunately, 
in such scenarios the spin-independent cross section for direct dark matter detection is predicted to be very 
small, cTp ^ < 10~^^ pb. Rates for collider production of the singlet-like hi will be very low. The x will appear 
in chain decays at the LHC and its roughly 100 GcV mass should be measurable with reasonable accuracy. 
However, to determine how singlet it is would require observation of a displaced vertex. Unfortunately, the 
predicted non-singlet content of the x for the SS scenarios is of order 1 — ~ few x 10~^ — 10~^, sufficient 
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to make the decays to the x prompt. 

Overall, it is clear that even a slight extension of the strongly constrained CNMSSM to allow non- 
universality for the Higgs soft-masses-squared and for opens up the phenomenological possibilities very 
considerably. One finds fairly good ideal-Higgs-like scenarios. In addition, the very intriguing singlet- 
scenarios that are consistent with all experimental constraints and give excellent flh^ become quite promi- 
nent. 
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